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ABSTRACT 



Context. Asteroseismology provides a unique tool for studying stellar interiors. Recently p modes have been detected on the bright solar-like 
star a Cen A thanks to high-precision radial-velocity measurements. However a better characterisation of these p modes is clearly needed to 
constrain theoretical models. 

Aims. We observed a Cen A during five nights using the HARPS spectrograph in order to improve our knowledge of the seismic properties of 
this star 

Methods. We performed high-precision radial-velocity sequences and computed the acoustic spectrum of a Cen A. 

Results. We identify 34 p modes with angular degree / = - 3 in the frequency range 1.8-2.9 mHz and amplitude range 13-48 cm s"', in 
agreement with previous seismic studies. We find an enhancement of the frequency scatter with the angular degree Z that indicates, considering 
the high inclination axis of a Cen A, rotational splitting and explains the low values of previously suggested mode lifetimes. We also derive 
new values for the small separations that take the effect of rotational splitting into account . 

Conclusions. Our seismic study of a Cen A leads to a list of now well identified p-mode frequencies and shows the importance of taking the 
rotation into account in order to properly characterise the p modes even in quite short campaigns. 

Key words, stars: individual: a Cen A - stars: oscillations - techniques: radial velocities 



1. Introduction 

The a Cen system is the closest stellar object to the Sun. It is a 
triple system with one M dwarf. Proximo Cen, in a large orbit 
around two stars closer to each other, a Cen A and B. Those are 
sun-like stars with spectral types, respectively, G2V and KIV. 

With V = -0.01 a Cen A is a very bright star Because 
it is part of a visual binary, good constraint s on its mass, 
M - 1.105 ± 0.007 Mq, have been obtained dPourbaix et al 



20021) . It is thus an ideal target for asteroseismology of sun- 
like stars. Several attempts to detect p-modes have beeri made 
dOelly et al."l986'; 'Pottasch et al.!ll992|; ISchou & Buz"asil200lh 
until Bouchy & Carrier (2002) identified twenty-eight oscil- 
lation frequencies using the CORALIE spectrograph. A five- 
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*Based on observations collected with the HARPS spectrograph 
at the 3.6-m telescope (La Silla Observatory, ESO, Chile: program 
075D-0800A) 



night multi-site run, coupling measurements from UVES and 
UCLES, led to the confirmation of thos e p-modes, as well 
as to the detection of additional ones (IButler et al.l 12004 : 
Bedding et alll2004 . 



Modelling a Cen A has proven challenging. It appears 
that stars with masses around 1.1 sta nd on the fringe 
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of ha ving a convective core -see e.g 
(120041) and references therein- which is very difficult 
to model properly. Thanks to the various available con- 
straints, oscillation 

and, hence, luminosity, effectiv e temp erature, metal- 
licitv (|Neu forge- Verheecke & Magain 1997), and radius 
(IKervella et al.l I2003h . se veral models have be e n computed 
using these constraints ( Thev enin et al.l |2002t ^T houl et aT 
2003t lEggenbergeret all 1200 4; MigHo & MontalbanI 12005 



Montalban & Mighd l2006h . However, it 



is difficult to find 
models that correctly account for all the observations simul- 
taneously as stated in these publications. Furthermore some 
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Table 1. Journal of radial velocity measurements. The min- 
imum and maximum S IN of spectra obtained each night is 
given. cr^MS represents the dispersion, expressed in m s 



Date 


Nb spectra 


Nb hours 


SIN 


o-RMi- [ms"'] 


2005/04/18 


835 


10.05 


150-550 


2.04 


2005/04/19 


1040 


10.59 


150-500 


2.10 


2005/04/20 


1038 


10.63 


250-500 


1.80 


2005/04/21 


992 


10.86 


100-500 


3.26 


2005/04/22 


1054 


10.69 


200-500 


2.44 



problems in the modelling of the interior of a Cen A are still 
partially unsolved, like the possible presence of a convective 
core. Therefore, better constraints on the p modes are needed. 

We observed a Cen A during a five night run with the 
HARPS spectrograph. This paper presents the results of this 
run. We first give an account of the observations (Sect. 2). In 
Sect. 3, we describe the techniques used to obtain the acoustic 
spectrum. We also explain the methods used for frequency ex- 
traction and identification. Thirty-four p modes were identified 
and for some of them we identified multiple frequencies. We 
discuss, in Sect. 4, the observed sc atter in the frequencies with 
respe ct to the asymptotic formula (' Vanda kurovll 1 967 ; Tassoull 
ll98Q) . We present evidence that this scatter shows the signature 
of rotational splitting, which contradicts a previous interpreta- 
tion involv ing pulsating modes with short lifetimes compared 
to the Sun (IKieldsen et al. 2005 ). New values are given for the 
small separations that take this effect into account. 



2. Observations and data reduction 

The target was observed over 5 nights in April 2005, using 
HARPS at the ESO 3.6 m telescope. We took sequences with 
typical exposures between 2 s and 10 s, with a dead time of 31 
s in between. In total, 4959 spectra were collected with typi- 
cal signal-to-noise ratios {S IN) per pixel in the range 300-450 
at 550 nm. The obtained spectra were extracted on-line and 
in real time using the HARPS pipeline. Wavelength calibra- 
tion was performed with ThAr spectra, and the radial velocities 
obtained by weighted cross-correlation with a numerical mask 
constructed from the Sun spectrum atlas. The journal of the ob- 
servations is given in Table 1 and the rad i al velo cities are pre- 
sented in Fig.[T] FoUowing lBouchy et al.l (1200 ih we computed 
the photon noise uncertainty on the radial velocity and found an 
average value of 18.5 cm s"' per measurement. The dispersion 
for each individual night is in the range 1.5-3.3 m s far from 
the photon noise value, and is dominated by acoustic modes 
with periods around 7 minutes, as shown in Fig. |2] Guiding 
errors, especially during nights 4 and 5, also introduce a non 
negligible source of noise as discussed in the next section. 



3. Acoustic spectrum analysis 

To compute the power spectrum of the velocity time series dis- 
played in Fig.[T] we used the Lomb-Scargle modified algorithm 




3 3 
BJD - S453479.5 



Fig. 1. Radial velocity measurements of a Cen A. The disper- 
sion for each individual night is given in Table [T] 
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Fig. 2. Zoom of radial velocity measurements in a three-hour 
sequence, taken during the third night, showing the presence of 
p-modes in the time series with periods of around 7 minutes. 
The semi-amplitude of about 3ms"' does not represent the in- 
dividual amplitude of p-modes but comes from the interference 
of several modes. 



( Lomb 19761 Scargle 1982 ) for unevenly spaced data. The re- 
sulting LS periodogram, shown in Fig. [3] exhibits a series of 
peaks between 1.6 and 2.9 mHz modulated by a broad enve- 
lope, which is the typical signature of solar-like oscillations. 
This signature also appears in the power spectrum of each in- 
dividual night. The power spectrum also exhibits three peaks at 
3.1,6.2, and 9.3 mHz. 

3.1. Guiding noise 

These peaks result from a recently identified guiding noise. 
This is caused by the telescope motor experiencing a hard point 
during its rotation, which has a period of approximately 5.5 mn. 
This problem affects the tracking system, causing the stellar 
light beam to move inside the optical fibre, without significant 
changes of flux. Such effects should be corrected by the guid- 
ing system. Unfortunately this can be done only if the exposure 
time is longer than the duration of one guiding cycle, approxi- 
mately 2 s, which was not the case during our run. Furthermore, 
the guiding camera exposure time, for very bright stars like 
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a Cen A, is as short as 0.01 second. With such a short expo- 
sure time and with good seeing conditions, speckles become 
visible and the guiding system no longer tracks the Gaussian 
centre of the stellar light spot seen by the fibre, but the bright- 
est speckle. Within these conditions the HARPS guiding sys- 
tem is not well-suited, so it provokes large displacements of 
the stellar spot inside the fibre. This behaviour of the guiding 
system is the reason for the radial velocity spikes we can see in 
Fig. [1] during the last two nights. They correspond to the best 
seeing conditions of the run. The field of view of the fibre is 1 
arcsec. In terms of radial velocity, the aperture of the fibre at 
the entrance of the spectrograph corresponds to a slit of width 
3 km.s"'. During commissioning tests, the typical scrambling 
was estimated around 400, which means that a spot moving 
from the centre to the edge of the fibre would induce a change 
in radial velocity of about 8 m.s"^ The magnitude of this effect 
decreases with the increasing size of the spot. With a 1 arcsec 
seeing, a 3-4 m.s"' effect on the radial velocity was measured 
during commissioning. In nominal cases, the RMS stability of 
the guiding system is about 0.1 arcsec, and the guiding noise 
is then estimated to 30-40 cm s"'. This noise is considerably 
reduced, and eventually vanishes, for long exposure times (>5 
mn), when the correction due to the guiding system cycle (2 s) 
is completely averaged. This was obviously not at all the case 
during our run. Especially during the two last nights, we clearly 
see several spikes with amplitudes from 4 to 8 m s"' in the RV 
times series. In the power spectrum the guiding system should 
introduce a significant contribution to the noise. 

Considering Gaussian noise, we found a 3.7 cm s ' mean 
noise level in the amplitude spectrum in the frequency range 4- 
5.5 mHz. Since the time series is based on 4959 measurements, 
the corresponding velocity accuracy is 1.47 m s"'. This value 
is far from the estimated photon noise of 18.5 cm s To verify 
this value, we separated the contributions to the radial velocity 
coming from the odd and even spectral orders of the spectro- 
graph. Because there is no dependence of the radial velocity on 
the parity of the order, subtracting the two curves will suppress 
all Doppler signals, including the contribution to noise due to 
the guiding. The remaining signal can be considered as pure 
photon noise. In the Fourier space we measure a noise of 0.5 1 
cm s indicating a photon noise of 20.3 cm s ' for the whole 
spectrum, in very good agreement with the first estimate. 

We can also compare this measurement to th ose issued 
from the HARPS run on n Ara (iBouchy et alJl2005h . The noise 
level was 1.17 m s which is significantly lower. But in that 
case the exposure time was 100s, allowing for the guiding noise 
to be averaged, at least partially. Consequently we may con- 
clude that the noise level of 1.47 m s"' comes in large part 
from the guiding noise and is seven times larger than the pho- 
ton noise. This analysis addresses the problems one can en- 
counter when observing bright stars with HARPS. It is neces- 
sary to use exposure times long enough to average the guid- 
ing contribution to the noise. In case this is not possible, as 
for a Cen A (because saturation would be reached with longer 
exposure times), it should be possible to reduce the effect by 
defocusing the telescope and hence increasing the size of the 
stellar light beam in the fibre. Another solution would consist 




Fig. 3. Amplitude spectrum of the radial velocity of a Cen A. 

in optimising HARPS guiding system with, for example, a tip- 
tilt system (.Bouchy & Connes,.1999,) . 

3.2. Echelle diagram 

In solar-like stars, p-mode oscillations are expected to pro- 
duce a characteristic comb-like structure in the power spec- 
trum with mode frequencies v„,; approxima ted reasonably 



trum witn moae rrequencies v„,; approxima tea reasonably 
well by the simplified asy mptotic relation ( Tassoull 1980t 
Christensen-Dalsgaardll988h : 



v„,/ ~ (« + - -H e)Avo - 1(1 + l)Do. 



(1) 



The two quantum numbers n and / correspond to the radial 
order and the angular degree of the modes. If the stellar disk is 
not resolved, only the modes with degree I < 3 can be observed. 
For given order and degree n and I, large and small separations 
are defined respectively by 

Av„,; = V„+i_; - V„j (2) 

and 

Sv„j = v„j - v„_i,/+2, (3) 
where Avq is the mean large separation, i.e. the mean frequency 
separation between two modes of consecutive radial orders, 
and is also defined by 

dr 



AvA 



Jo 



(4) 



with R the stellar radius and c the sound speed. Do is given by 

< Sv„j >„ 



Do 



(5) 



4/ + 6 

Frequencies separations, as well as frequencies of ind i vidua l 
mo des, were previ ously reported by iBouchv & Carried ( 120021) 
and lBedding et alJ §004). 

The power spectrum shown in Fig. [3] presents a clear and 
unambiguous periodicity of 106 juHz both on its autocorrela- 
tion or in the comb response, giving an estimate for Avq. To 
identify the angular degree / of each mode individually, we di- 
vided the power spectrum into slices of 106 juHz and summed 
them up. The corresponding summed echelle diagram is pre- 
sented in Fig.|4]and allows us to unambiguously identify modes 
/ = 0, 1 , 2 and their side-lobes due to the daily aliases at + 1 1 .57 
juHz. Some low-amplitude / - 3 modes were also identified. 
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20 40 60 80 100 

Frequency modulo (106 /^Hz) [/xHz] 



Fig. 4. Sum of the echelle diagram of a Cen A. Peaks corre- 
sponding to the / - 0,1,2 and 3 and their daily aliases are 
identified (with a prime). 



16 
17 
18 
19 

20 
21 
22 
23 
24 
25 
26 



: 



/= 1 



1 = 1 



1837.6 
1943.0 

2152.6 

2363.3 
2470.4 
2576.3 
2682.^ 
2788.zl 
2895.fi 



1886.5^ 
1990." 

2096.8 

2202.0 

2308.7 

2415.^ 

2519.6 

2626.4 

2840.1 



1936.6^ 
2038." 

2143.9^ 
2250.8/2254.4 
2355.8/2357.6 
2462.6/2465.0 
2568.^2571.8 
2674.4/2676.S 



'Only with gap-filling 

^Only on-sight (with or without gap filling) 

*Only with CLEAN 



/ = 3 



2299.4 



2511.5 



2833.8 



3.3. Identified frequencies 

We distinguish two different stages in the mode identification 
process. First we need to extract the frequency from the acous- 
tic spectrum. Roughly speaking this phase consists in locating 
the high-amplitude peaks in the spectrum. Second we have to 
verify whether it indeed coiTesponds to a pulsation mode of 
the star We used two different methods to extract the frequen- 
cies from the amplitude spectrum, as described below. As a 
general rule, among the peaks extracted from the spectrum, we 
selected only the strongest, i.e. with S IN ratios higher than 3 
(> 11 cm s '). However, we needed some additional selection 
process to determine whether each frequency was in acceptable 
agreement with the asymptotic relation. This kind of process is 
characteristic of frequency identification in Sun-like stars. It is 
equivalent to assigning a true or false value to the peaks we can 
extract from the amplitude spectrum. 

It is worth noting that, regardless of the extraction method 
used, the spectrum can be treated to improve the mode detec- 
tion. The coupling between the gap-filling technique and the 
on-sight extraction, as described below, is an example of such 
an additional process. 



Table 3. Amplitudes of the identified p-modes (cm s '). 



We first used a CLEAN algorithm (Robe rts et al.lll987h . 



This method presents the advantage of removing the artifacts 
in the spectrum resulting from the gaps in the time series. 
However, after using CLEAN, we still needed to reject some 
frequencies that were obviously false p-modes identifications. 
This raised the question of the limits of the CLEAN procedure. 
The main problem concerns the finite lifetimes of the modes: 
using CLEAN, we subtract pure infinite sine waves from the 
signal, despite the fact that the observed modes have a finite 
lifetime. There is a risk of removing signal that was not origi- 
nally in the spectrum. 

The second method supposes an on-sight extraction of the 
frequencies in the amplitude spectrum. We select each fre- 
quency with respect to the template obtained from the summed 
echelle diagram (Fig. Hjl. Unlike the CLEAN algorithm, it 
does not modify the original spectrum. We also note that 
the extraction and frequency selection are simultaneous. The 
limits of this method are obvious since many structures in 





/ = 


/= 1 


1 = 2 


/ = 3 


n= 16 


21 


13 


13 




n = ll 


14 


24 


22 




n = 18 




32 


12 




n=l9 


22 


15 


52/26 


16 


n = 20 




37 


13/30 




n = 2l 


20 


22 


22/23 


14 


n = 22 


19 


48 


25/20 




n = 23 


23 


25 


29/18 




n = 24 


16 






16 


n = 25 


13 


14 






n = 26 


14 









the spectrum cannot be distinguished on sight and it cannot 
account for interferences that may occur between different 
modes. We also proceeded to an on-sight identific ation using a 



pre-treated spectrum with the g ap-filling method (iFossat et al 



19991: iBouchy &Carrieill2002l) . This method aims at reduc- 
ing the efi'ect of single-site observation. It consists in replacing 
gaps in the time series by signals obtained at a time t = 2/Avo, 
before or after It is justified since a peak appears in the auto- 
correlation spectra for Avo/2 (.Bouchv & Carrier,200Z) . 

In all cases, n values are assigned to p modes using the 
asymptotic relation ([T]i. We assumed that the parameter e has a 
value close to the solar one (e© ~ 1.5). The frequency and am- 
plitude of the 34 p-modes we identified are listed in Tables |2] 
and[3] Generally the two extraction methods are in good agree- 
ment. The idea behind using them was to check their agreement 
and to obtain complementary sets of frequencies. Therefore, 
in the following work we used all the frequencies shown in 
Table|2]and we did not assign weights depending on the detec- 
tion technique used. However, the frequencies not appearing 
for all the techniques are marked in Table |2l For some I = 2 
modes, we identified multiple frequencies. 

Amplitudes were determined by assuming that none of the 
p modes are resolved and that it corresponds to the height of 
the peak in the power spectrum after quadratic subtraction of 
the mean noise level. Considering that the frequency resolution 
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30 40 50 60 70 
P'rcqucncy modulo 106 j.i\Az 



Fig. 5. Identified p-modes /=0 (o), /=! (.), /=2 (a), and Z=3 (A) 
in echelle diagram format. The size of the symbols are propor- 
tional to the amplitude of the modes. The crosses correspond 
to ambiguous modes (see text for details). 



First the asymptotic relation in the form given by Eq. ([TJ 
is valid for high fr equencies, i.e. large compared to t he Brunt- 
Vaisala frequency ( Vandakurovll 1 967 : Tassoul Il980h . There is 
a small departure from this relation in the echelle diagram. 
However, for the considered radial orders in the range 16- 
26, this approximation is accurate. This can be seen in Fig. [8] 
where the large separation is nearly constant with increasing 
frequency. 

Second, other physical processes can add some scatter 
around each individual frequency, v„j. In the first instance, be- 
cause the oscillation modes in sun-like stars are stochastically 
excited and have finite lifetimes, they appear in the power spec- 
trum as a series of peaks modulated by a Lorentzian. The half 
width of the Lorentzian reflects the dampi ng rate of the mode 
and thus its lifetime, see e.g. lToutain & Froehlich (I1992I) . This 
effect has been well-studied for the Sun, leading to characteris- 
tic lifetimes betw een 38.4 and 3.15 day s in the frequency range 



1570-2808 ^Mz (IChaolinet al.Ml997h . iKieldsen etal.1 (120051) 



suggest that the scatter they measured around the frequencies 
in a Cen A correspond to this finite lifetime effect. They give 

days at 2.6 



the values of 2. 
mHz. 



+0.9 



yj^l days at 2.1 mHz and 2.1+;j 



Furthermore, since the star is rotating, the proper frequency 
of a given mode will be split, in the observer frame, following 
the relation 



of our time series is 2.6 /iHz, we adopted an uncertainty on no- 
resolved oscillation modes of 1.3 //Hz. Such an uncertainty can 
be considered as conservative in the case of infinite time-life 
modes with high S /N ratio (>3.5). 

All the identified modes are displayed in the echelle dia- 
gram in Fig.|5] The inferred values for the radial order, n, and 
the angular degree, I, are also given in Tables |2] and [3] In some 
cases, we observed two frequencies corresponding to the same 
couple (n, 0- These multiple frequencies are given for I = 2. 
We note that, for each multiple detection at / = 2, at least 
one frequency was detected simultaneously with CLEAN and 
on sight. It therefore indicates that these frequencies are not 
just the same one identified slightly differently when changing 
the extraction method. We also extracted peaks at frequencies 
1987.0 and 2412.8 juHz. The first one cannot be identified pre- 
cisely as an Z = 1 or Z = 3 mode. The second one is a pecuUar 
case. If we select it, we obtain a multiple detection for an Z = 1 
mode. However, the situation is quite different from our / - 2 
cases, since both components here have been identified using 
a different method. It is thus possible that these two values re- 
flect a difference in the extraction methods used. Therefore, we 
decided not to include these two frequencies in Table|2l but we 
still display them in Fig.|5] 

3.4. Frequency scatter analysis 

If the frequencies were to follow the simplified asymptotic rela- 
tion given by Eq. ([1]) exactly, their representation in the echelle 
diagram would consist of parallel vertical lines, one for each 
degree. Although the observations reproduce this general trend, 
we observe a more complex structure in Fig.|5] 



V«,/,m - VnJ + mQ. 



(6) 



where m is the azimuthal order (-/ ^ m ^ I) and Q. the rotation 
rate. We suppose a rigid b ody rotation here. In the case of the 
SunJ Chaplin et al. ( 200lh find an average value of 0.4 //Hz for 
the rotational splitting at low degre es (/ < 3). Following their 



detection of p-modes in or Cen A, Bouchy & Carried (I2OO: 
suggest that they could have observed hints of a rotational split- 
ting of the order on 0.5 //Hz, compatible with an almost solar 
rotational period. 

Scatter analysis is a very interesting tool, especially when 
dealing with short seismic runs that do not provide a frequency 
resolution high enough to allow more direct t echniques, such as 
direct fitting of the acou s tic power spectrum (IGizon & Solanki 
2OO3I: iBallot et alj|2006t iFletcher et alj|2006l) . To measure the 
scatter, we first needed to estimate the departure from the 
asymptotic relation for each frequency, since we used it as a 
reference for the mode identification. For this purpose, we fit- 
ted the frequencies, for a given degree, with seco nd-order poly 



nomia ls. This method is close to the one used by lKjeldsen et al 



(l2005h . However, since we observe multiple frequencies (and 
do not average them) for some given (n, I) values, we could 
not use the same variables as they did (see their Eqs. (2)-(5)). 
We instead adopted the variables displayed in the echelle dia- 
gram, i.e. the frequency and the frequency modulo the average 
large separation, the latter one being a second-order polynomial 
function of the former one. 

The scatter cr*'' on the frequencies for a given degree /, with 
respect to the corresponding polynomial fit, is defined as the 
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10 20 30 40 50 60 70 80 90 100 
Frequency modulo 106 /v-H/ [/^-H/] 

Fig. 6. Echelle diagram for all the identified frequencies in 
a Cen A provided with HARPS (black-square), CORALIE 
(open-square) and UVES-UCLES (open-circle). The lines cor- 
respond for each angular degree to the fit of a second-order 
polynomial to all frequencies. We used them as guidelines to 
compare the different sets of frequency. 

square root of the normalised weighted sum of the fit residuals: 

( N, ' 



(M-3)^fl,- 



!=1 

with f"'°if/l) the rth frequency observed at degree / modulo Avo, 
a, the corresponding amplitude, f"y/^(l) the fitted value, and Ni 
the total number of frequencies observed at /. As a first approxi- 
mation, we consider that the amplitudes of the modes represent 
suitable weights. 

Assuming normally distributed frequencies around the fit- 
ting polynomial, the uncertainties on the computed scatter can 
be approximated by cr^'* / VZ^V, with the number of points 
used for the computation. The resulting scatters for individual 
degrees ar4i|cr|J^p = 0.50 + 0.12 //Hz, cr^^ = 0.70 + 0.17 //Hz, 
ando-[J^p = 1.71 + 0.34 //Hz. 

The I - 3 modes were excluded from our analysis. We only 
observed three of them and, thus, cannot derive a proper scatter 
using a second-order polynomial fit. All these modes have low 
amplitudes, so discarding them should not alter our analysis. 

As discussed in Sect. 13.51 this scatter increase with mode 
degree can be a consequence of rotational splitting. We 
searched for such trends in the previously published sets of fre- 
quencies. In Fig. |6] we represent the values from the previous 



*In the following, the subscript HRP will stand for HARPS, COR 
for CORALIE, and UVC for UVESAJCLES 



papers in the echelle diagram, in addition to our frequencies. 
The lines appearing in this echelle diagram are fits to all the 
frequencies and we used them as visual guidelines for the com- 
parison between frequency sets. It is worth stressing that the 
general agreement between these three sets is very good. 

CORALIE frequencies. From the pubHshed frequencies, the 



scatters are cr, 

J2) 



.(0) 
COR 



0.96 ± 0.24 //Hz, cr] 



,(1) 

COR 



0.90 + 0.20 //Hz, 



and cr^ojj - 1.01 + 0.23 //Hz. These show a marginal increase 
with degree /. However, comparing the three sets of frequencies 
in Fig . |6] suggests that the identifications in iBouchv & Carrier 
(I2OO2I) for the Z = modes at frequencies 2573.1 //Hz, 2679.8 



//Hz, and 2786.2 //Hz should be modified. We assign them a 



degree / = 2. The resulting scatters are o"cor - 0.33 ±0. 11 //Hz, 



COR 



0.90 + 0.20 //Hz, and cr 



,(2) 
COR 



1.53 + 0.30 //Hz. In 



this case the trend towards increasing at higher degrees clearly 
appears. This is not surprising since we obtain two possible 
new multiple frequencies for Z = 2, at radial orders 22 and 24, 
when changing the identification. 

UVES/UCLES frequencies. The picture for the 
UVES/UCLES frequency set is less clear. It should first 
be noted that multiple frequencies were also found for certain 
(n.l) values in iBedding et al.l (12004 ,). Nevertheless, in order to 
retain only one frequency, those multiple values were averaged 
using a weighted-mean. We also note that, based on their ridge 
identification technique, they retained modes with S /N <3. In 
order to compare homogeneous sets of frequencies, we only 
retain their modes with S /N > 3 (i.e. with amplitudes larger 
than 8.7 cm s"'). As a first approximation we can consider that 
the effect of averaging the multiple frequencies is negligible. 
To support this assumption, we raise the fact that all the mul- 
tiple frequencies identified are composed of a high-amplitude 
and a low-amplitude peak. The weighted mean is dominated by 
the latter one. Furthermore, the mode at 2572.7 //Hz has been 
identified as both I = and Z = 2 in their Table 1 . Comparing 
with the HARPS and CORALIE frequencies in Fig. |6l we 
choose to retain only the second identification. The resulting 



scatters are cr|°' = 1 .04 + 0.25 //Hz, (t\1> = 1.32 + 0.28 //Hz, 



and cr. 



(2) 



UVC 
1.32: 



,(1) 

UVC 



' UVC ~ i •-'^ 0.28 //Hz. The trend is much more marginal 
in this case and only visible between I - and / = 1 scatters. 
The fact that multiple frequencies are observed for both / = 1 
and / - 2 modes may be the main reason for this lack of 
increase. This may be coupled to resolution effects (see discus- 
sion in Sect. 13. 5t . The same kind of behaviour is observed if 
we compute the scatters from the non-averaged frequencies (T. 



Bedding, private communication), cr. 



UVC 



1.84 ± 0.38 //Hz, and cr] 



(2) 
UVC 



.(0) 
UVC 



1.13 ± 0.27 //Hz, 
1.77 ± 0.36 //Hz. We 



should eventually point out that if we compute the scatter while 
retaining all the published frequencies (i.e. without cutting at 
S /N > 3), it decreases between / - 1 and I = 2 modes. 

The above discussion shows that, for the HARPS and 
CORALIE frequency sets, the scatter shows a clear trend to- 
wards increasing with the mode degree. It is not possible to 
distinguish such a trend in the UVES/UCLES frequencies, al- 
though the large uncertainties on the scatters for / - 1 and 
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I - 2 do not exclude it. It is also worth stressing the sensi- 
tivity of the scatter determination to the mode identification, 
as shown with the CORALIE and UVES/UCLES frequencies. 
This observed trend could nevertheless help for characterising 
the multiple frequencies observed for I - 2 modes. 

3.5. Characterisation of the split modes 

The m-dependent visibility of a rotationally split mode (I, m), 
independently of the radial order «, is linked to the incli- 
nation ; of the rotatio n axis of the star by the coefficient 
dOizon & Solankill2003b : 



(i-\m\y. ^ 

(l + \m\)l I'' 



[pl'"l(cosO]^ 



(7) 



where f is the Legendre polynomial. 

According to Eq. (|7]) at high inclinations (i.e. / > 65°), 
modes with the highest visibility coefficients <?/„,(/), and thus 
the highest amplitudes, are those with / = \,m - 1 and / = 2, 
m - 2. Modes with I - 2, m - Q could possibly be seen, but 
their amplitude in the power spectrum is expected to be approx- 
imately two times lower than their m - 2 counterparts. Modes 
with I - \,m - Q and / = 2, m = 1 almost disappear for high 
inclinations. 

In the framework of the rotational splitting hypothesis, we 
can interpret the multiplets found in Table |2] for some / - 2 
modes as rotationally induced. This assumption is justified by 
the above discussion. If we make the assumption that the incli- 
nation axis of a Cen A is aligned with the axis of the binary 
system, we can then consider ; - 79° ( Pourbaix et alj [l999). 
It follows that modes with (/, m) - (I, I) and (2,2) should be 
the most visible. Whereas the effect of mode lifetimes on the 
scatter should not depend on the mode degree, an enhancement 
with increasing I and, therefore, \m\ is possible, explaining the 
trend observed in HARPS and CORALIE scatters. We eventu- 
ally note that assigning azimuthal orders to the split modes we 
observe does not exclude a contribution coming from the finite 
mode lifetime effect to the overall scatter, as discussed in the 
following. 

The measured separations in our frequency multiplets at 
/ = 2 are plotted in Fig.|2l The error bars have been calculated 
using the quadratic sum of the uncertainties on our frequen- 
cies. The largest possible splitting effect would be obtained 
for modes with m - +2. Assuming an almost solar rotation 
rate, the difference would be around 2 /^Hz, which is com- 
patible with the values displayed in Fig. [T] However, the di- 
agnostic potential of these separations should be considered 
cautiously. On one hand, if the mode lifetimes are long com- 
pared to the observing time, T, then the accuracy to which two 
close frequencies can be resolved is limited by our resolution, 
which scales approximately like 1 jT . Depending on the phase 
difference between the modes, they could result in only one 
peak in the power spectra or, on the contrary, two peaks sep- 
arated from > I IT, even if the true distance between them 
is smaller than the measured splitting (jLoumos & Deeming 



I978l:IChristensen-Dalsgaard & Goughlll982l) . This could lead 
to overestimate the rotational splitting. We note that, under 



3500 
Frequency [/xHz] 



Fig. 7. Measured splitting for / = 2 modes in HARPS frequen- 
cies. Symbols are the same as in Fig.|5] 



the above hypothesis on the mode lifetimes, this resolution ef- 
fect could be involved in the large scatter we observed in the 
UVES/UCLES frequencies at / - \. This effect certainly im- 
plies that a small distance in frequency in the power spectrum, 
caused by splitting for / - 1 modes, will be more severely bi- 
ased than a larger one, caused in this case by an / = 2 splitting. 
On the other hand, if the mode lifetimes are similar to or shorter 
than the observing time, their contribution to the scatter could 
become similar to or greater than the effect of rotation. In this 
case, it becomes difficult to disentangle the two contributions 
and heavier simulations are needed. 



Gizon & Solankil ( 120031) and lBallot etal.1 (12006b used sim 



ulated seismic data, spectra corresponding to long time series 
of respectively 180 and 150 days, to study the possibility of es- 
timating the rotation rate and the inclination of a sun-like star. 
The spectra models used to describe these simulations were as- 
signed Q. and / as the only two free parameter s. These were esti- 
mated using a maximum-likelihood method. iGizon & Solanki 
( 2003) cons i der a single-split mode in the spectrum, whereas 



Ballot et al.l (12006) use a multi-mode fitting approach. In rela- 
tion to this study, the main conclusion from both studies is that, 
for a rotation rate Q = Oq, it is not possible to determine ac- 
curately the inclination that introduces a bias in the estimated 
rotational splitting. Moreover, because of computational costs, 
their studies were made with fixed mode lifetimes. It is not 
clear which degree of complexity would be added to this prob- 
lem if mode-lifetime estimation were to be included. Given the 
length of the time series considered in these works, it is doubt- 
ful that a conclusion could be reached using short seismic runs 
such as the one presented here. Still, this gives a good insight 
into the overall problem of estimating stellar parameters from 
acoustic spectra of solar-like stars. It also points out the impor- 
tance of having an independent estimate of the inclination, as 
in our case, in order to estimate a rotation rate using asteroseis- 
mology. 



Kjeldsen et al.l (l2005h estimated mode lifetimes in a Cen A 



using the UVES/UCLES frequencies. Their method consists 
in simulating several realizations of a time series correspond- 
ing to their observations and from them extracting frequencies 
to obtain a theoretical estimate of the scatter. They needed to 
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compute at least 100 simulations per set of input parameters. 
However, they supposed a static star in their study, and thus 
neglected rotation. Using this technique with an additional pa- 
rameter to estimate is possible, although it would represent sev- 
eral months of CPU time. Describing possible methods to by- 
pass this prohibitive computing time is beyond the scope of 
this paper and will be addressed in a later study. Another con- 
sequence of our identification of splitting is that neglecting ro- 
tation may lead to underestimating mode lifetimes, since they 
should not account for all th e observed s catter. The value given 
for the mode lifetimes in iKieldsen et all fcoOSi) should thus be 
considered as a lower limit. 

Independent measurements of the rotation of a Cen A 
have been made, which could indicate the expected splitting. 
A rotational period of 28.8 + 2.5 days has been measured 
for a Cen A from analysis of the modulation of lUE spectra 
dHallam et al ][l99ll) . However, this paper states that several so- 
lutions come out of the analysis. Although the result they give 
is their best estimate, they are not able to definitely rule out 
other shorter rotational periods. The splitting co rresponding to 
this ro tational period would be 0.40 + 0.04 juHz. ISaar & Osten 
( Il997h derived a projected velocity, vsin/ - 2.7 + 0.7 km.s~\ 
from high-resolution observations of Fel, Cal and Nil lines. 
Using the value of the radius determined by interferometry, 
R = 1.22 4 ± 0.003 Rrr. tKervellae t al. 2003) and the incUna- 
tion from lPourbaix et alJ (1 19991) . we obtained P = 22.5 + 5.9 
days from this measurement and a corresponding splitting Q. - 
0.51 + 0.13 juHz. 

More recently. lFletcher et al. I (l2006l) used a 50-day time se- 
ries from the WIRE satellite to estimate the rotational splitting 
and mode lifetimes of a Cen A. They used both power spec- 
trum and autocovariance fitting techniques. Their best estimate 
gives a rotational splitting of 0.54 ± 0.22 yuHz and a mode life- 
ti me of 3.9 ± 1 .4 d a ys. Th is value is higher than the one derived 
by iKjeldsen et alJ (l2005h . which supports our analysis stating 
that the mode lifetime should not be the sole contributor to the 
frequency scatter and that neglecting the splitting should intro- 
duce a significant bias. 

3.6. Large and small separations 

Individual large and small separations obtained from the new 
HARPS frequencies are represented in Fig.|8]and|9]and uncer- 
tainties on both of them are 1.8 yuHz. For / = 2 modes with 
multiple frequencies detected, we used the mean frequency 
value in Eqs. (|2]l and (|3]l. The small separations are very use- 
ful for stellar modelling and ne ed to be strongly constrained, 
e.g. lMontalban & Miriiol ( l2006l) . The first two HARPS small- 
separation values were derived using I = 2 modes with only 
one frequency. It should be emphasised that, considering rota- 
tional spUtting, these values are questionable since one m = ±2 
mode could be missing. 

Fig.|9]also displays the small separations obtained from the 
CORALIE and UVESAJCLES runs. We note that the HARPS 
values do not show the trend of CORALIE small separations 
towards decreasing sharply at high frequencies. In that re- 
spect, they show very good agreement with the separations de- 



3500 
Frequency [/j-Hz] 



Fig. 8. Large separations Av obtained from HARPS frequencies 
versus frequency for p-modes of degree l-Q, 1 and 2. Symbols 
are the same as in Fig.|5] 



2000 3500 
Frequency [/j-Hz] 

Fig. 9. Small separations 6v„i for I = 0. The dififerent symbols 
correspond to the values obtained with HARPS (black-square), 
CORALIE (open-square) and UVES-UCLES (open-circle). 
Crosses mark the two CORALIE small separations which dis- 
appear after the new identification suggested in Sect. 13.41 For 
sake of clarity the error bars were omitted. 



rived from the UVESAJCLES run. Bedding et al.1 (120041) also 
identify some modes with multiple frequencies. They com- 
pute the corresponding separations using a weighted mean 
(using the amplitudes as weights). Although their treatment 
diff'ers slightly from ours, we consider it a small eff'ect that 
does not challenge the overall agreement between HARPS and 
UVES/UCLES small separations. Eventually, the new identi- 
fication suggested in Sect. 13.41 for the CORALIE frequencies 
results in the disappearance of the last two small separations 
at 2573.1 yuHz and 2786.2 //Hz for this set of frequencies. We 
note that this improves the agreement between the three runs. 

We note that averaging without weighting our I = 2 fre- 
quencies may lead to underestimate the contribution to the fre- 
quency scatter coming from the mode lifetime. On the other 
hand a weighted average is unnecessary for rotationally-split 
modes since their contribution to the mean frequency should 
be the same when averaging over m. The fact that the HARPS 
and UVES/UCLES small separations agree well shows that the 
assumption made before averaging does not have a huge effect 
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on the computed values. The important point is therefore to ob- 
tain more than one frequency, at least for I - 2, so that both 
effects can be taken into account by using mean frequency val- 
ues. 

Furthermore, from the results displayed in Figs.[8]and|9] we 
obtained new mean values for the large and small separations 
for the HARPS frequencies, respectively, 105.9 ± 0.3 yuHz and 
6.9 + 0.4 yuHz. The latter value was obtained by averaging the 
small separations over n. The former was estimated by linearly 
fitting the relation v„j - nAvo + k, because consecutive values 
of the large separation are not independent; here ^ is a constant 
and the slope is the average large separation. Using the same 
methods, we obtained 105.5 ± 0.3 yuHz and 6.8 ± 0.4 yuHz for 
CORALIE (using the new identifications given in this paper) 
and 106.1 ± 0.4 //Hz and 7.1 ± 0.6 //Hz for UVES/UCLES. 
These values stress the very good agreement existing between 
the different runs. 

4. Conclusion 

We analysed observations of a Cen A made during a five- 
night run using the high-precision spectrograph HARPS. We 
extracted 34 frequencies from the amplitude spectrum, con- 
firming the previous studies of this star and were able to im- 
prove the characterisation of the observed p-modes. 

Our main result came from analysing the scatter in these 
frequencies. We found clear evidence of rotational splitting in 
the spectrum of a Cen A: an increase in the scatter with the 
angular degree of the modes and multiplets for I - 2 modes. 
This observation suggests that the valu e of the mode Ufetime 
should be higher than the one given by iKjeldsen et al.l (l2005h 
without taking the rotation into account. We thus conclude that 
a frequency scatter analysis must include the rotational splitting 
effect. 

Observing rotational splitting allows us to derive precise 
values for the large and small separations, using frequen- 
cies averaged over the azimuthal order m. These new val- 
ues are in good agreement with t he small separations ob- 
tained from the UVES/UCLES run (IBedding et al.ll2004 !) and 
contradict the trend seen in those obtained from CORALIE 



(iBouchy & Carried l2002h . This result is fundamental for fur- 
ther modelling of a Cen A. 

Although we argue against short mode Ufetimes, this prob- 
lem still requires more observing time, on longer periods, to be 
fully addressed. Such observations would also allow us to ob- 
tain a better resolution, which is necessary for detecting very 
close frequencies unambiguously and obtaining direct mea- 
surements of the rotational splitting. 

This study shows that ground-based, single-site, seismic 
observations are valuable for improving our knowledge of the 
pulsation modes in sun-like stars. Stronger seismic constraints 
are also needed on a Cen B in order to model the entire sys- 
tem efficiently. These results are also very encouraging in the 
framework of future base-ground experiments. It is highly rel- 
evant that rotational splitting could be measured even from 
single-site measurements and demonstrates the strong possi- 
bilities offered by asteroseismology, especially if longer times 
become available. 
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